Photocatalytic-ozonation process in oxytetracycline degradation in aqueous solution: composite characterization, optimization, energy consumption, and by-products

In this research, we synthesized BiOI/NH2-MIL125(Ti) via solvo-thermal method to investigation of oxytetracycline (OTC) degradation in photocatalytic-ozonation process. The results of the XRD, FESEM, EDAX, FTIR, UV–Vis, TEM, XPS, and BET analyzes indicated that the catalyst BiOI/MOF was synthesized with excellent quality. Design of experiment (DOE), ANOVA statistical analysis, interaction of parameters and predicated optimum condition was done based on CCD. The effect of catalyst dose (0.25–0.5 mg/l), pH (4–8), reaction time (30–60 min) and O3 concentration (20–40 mN) at 10 mg/l of OTC on PCO/O3 process was optimized. Based on P-value and F-value coefficients (0.0001, 450.3 respectively) the model of OTC (F-value = 2451.04) and (P-value = 0.0001) coefficients, the model of COD removal was quadratic model. Under optimum condition pH 8.0, CD = 0.34 mg/l, RT = 56 min and O3 concentration = 28.7 mN, 96.2 and 77.2% of OTC and COD removed, respectively. The reduction of TOC was 64.2% in optimal conditions, which is less than the reduction of COD and OTC. The kinetics of reaction followed pseudo-first-order kinetic (R2 = 0.99). Synergistic effect coefficient was 1.31 that indicated ozonation, presence of catalyst and photolysis had a synergistic effect on OTC removal. The stability and reusability of the catalyst in six consecutive operating steps was acceptable and 7% efficiency decreased only. Cations (Mg2+, and Ca2+), SO42− had no influence on performing the process, but other anions, organic scavengers, and nitrogen gas, had an inhibitory effect. Finally, the OTC degradation probably pathway includes direct and indirect oxidation that decarboxylation, hydroxylation, demethylation and were the main mechanism in OTC degradation.

www.nature.com/scientificreports/ molecules, and semiconductor catalysts in combination. Due to its attractive properties such as inexpensive, eco-friendliness, availability, suitable efficiency, and simple operation compared to other processes, it creates a fundamental function in the organic matter degradation 10 . It was shown that this integrated process enhances reactive oxygen species and subsequently enhances mineralization rates of resistant organic compounds, such as OTCs. The ozone molecule has just 2.07 V standard oxidation-reduction potential (ORP), while PCO/O 3 process uses free oxidation radicals such as hydroxyl radical and reactive oxygen species (ROSs), that have higher than 2.0 eV ORP. In the photocatalytic process (PCO), the OTCs decomposition is restricted close the catalyst's surface, whilst in the integrated process, OTCs can be degraded to inorganic compounds such CO 2 and H 2 O via ozone in solution 11 . The PCO/O 3 process mechanism described by (Eqs. [1][2][3][4][5][6][7][8][9][10][11][12] 10 . At the begging of the reaction chain (Eq. 1), UV light irradiates cause of the production of electrons and holes. Based on (Eq. 2) reaction, ozone molecule uses of photo-generated electron to generate O 3 ·− . The oxygen captures electrons form HO 2 · /O 2 · , that could react with ozone and generate O 3 ·− . According to (Eq. 3), O 3 ·− reacts via H + to be transformed into HO 3 ·− and OH · (Eq. 4). Based on (Eq. 5), Hydrogen peroxide (H 2 O 2 ) generated that produces OH · and O 3 ·− based on subsequent reactions. Newly, study about PCO/O 3 process focalized in presentation of novel semiconductors. Some semiconductors like metal-organic framework (MOFs), g-C 3 N 4 , and slim band-gap oxide of metals families were used to expand the response spectrum of light until visible range 12 . The classic semiconductors like TiO 2 and ZnO activated via UV-C (under 300 nm wavelength that is only 4% of sunlight spectrum) 13 . Recently, a growing attraction in MOFs has focused on photocatalysis for the photo-driven removal of pollutants 14 . There have been a variety of applications of them in recent studies, such as catalysis, chemical sensing, and adsorption 15 . Various MOF-based photocatalysis have recently been developed in environmental studies. The Ti-based amino-functioned MOF (NH 2 -MIL125 (Ti)) was chosen for some advantages, including non-toxicity, being cheap, photo/water stability, and visible-light absorption 16 . The BiOI/NH 2 -MIL125(Ti) catalyst was prepared to make a synergistic effect with outstanding photocatalytic activity 17 . Bismuth-based semiconductors as a novel form of nonpoisonous, consistent photocatalysis via high visible light absorbency are worthwhile, environment-friendly, and stable compounds 18 . Among these, BiOI has shown impressive photocatalytic activity due to its appropriate band gap (BG = 1.62-1.93 eV), exclusive layered structure, and Special electrical features. The high rate of recombination of BiOI restricts its practical application. It has been proposed that heterojunctions can overcome this disadvantage by increasing charge separation and reducing the photo-generated electrons and hole recombination 19 . One way to improve the efficiency of photocatalysis is to create a heterojunction structure. This method reduces the energy of the band-gap between the conduction band and the valence band, and as a result, the photocatalyst is activated in the range of visible and UV-A light. Also, the structure of catalyst BiOI/NH 2 -MIL125(Ti) is such that it overcomes the limitations of structures BiOI and NH 2 -MIL125(Ti) and finally becomes a catalyst with proper efficiency in the PCO/O 3 process 20 . Alexsandra Valério et al., an investigation of the photocatalysis and ozonation process on the tetracycline degradation. In this study, the total organic carbon (TOC) removal reached above 90% after 3 h. Compared to the sum of the separated processes, the kinetics of reaction increase by 20% 21 . In this current study, for the first time, BiOI/NH 2 -MIL125(Ti) was used as a catalyst in the PCO/O 3 process for the decomposition of OTC. This catalyst possesses a remarkable capacity to generate active free radicals because of its unique (12) Free radicals + OTCs → CO 2 + H 2 O + non -toxic inorganic compounds where D represents the crystal samples size, λ is the angle of diffraction, K represents the dimensionless constant, β is FWHM (full width half maximum) of the diffraction peak and θ is the diffraction angle and θ is the diffraction angle.
The optic properties and structural features were examined using UV-Visible spectrum (UV-Vis DRS) was recorded by Agilent Cary 60 spectrophotometer. Tauc equations (Eqs. 14, 15) were used to determination of band-gap energy 25 .
where, α is the coefficient of absorption, v is frequency of light, A is a proportionality constant, h is Planck's constant, and E g is the band-gap. The band-gap energy (direct band-gap semiconductor) of BiOI is calculated using Eq. (14) while that of MOF (indirect band-gap semiconductor) is calculated using Eq. (15). The morphology of samples was observed by FE-SEM (UN41219SEM) under vacuum conditions of ≥ 1.2 × 10 −4 mbar. An analysis of elemental mapping and purity was performed on samples using energy dispersive spectroscopy (EDS). Shimadzu JEM-1200 EX with 100 kV stimulating voltage was used to obtain transmission electron microscopy image (TEM). Nitrogen adsorption at 77 K was used to determine the volume, pore size, and surface area distribution of BiOI/MOF. In situ degassing of samples was conducted for 12 h at 200 °C under vacuum. Brunauer-Emmet-Teller (BET) applied to calculate surface areas according to the linear relationship between p/ p 0 and surface areas. An ESCALAB250XI electron spectrometer was used to analyze the XPS characterizations. www.nature.com/scientificreports/ Photocatalytic ozonation process set up. PCO/O 3 process runs were done in a 600-ml cylindrical unit with D = 10 cm and h = 20 cm dimensions. The quartz sheath with dimensions 4.5 × 16 cm was installed horizontally in the middle of the reactor. The UPVC box with 39.5 × 41.5 × 30 cm was used as a pilot unit (Fig. 1).
A batch system was used to perform the runs at ambient temperature (25 ± 2 °C). 0.1 g of OTC was dissolved into deionized-water for preparation of stock solution (100 mg/l). The UV-A (Low-pressure, λ = 385 nm maximum wavelength range, and 6 W) was installed in the middle of reactor, with 15 cm length (TUV-PLL6W-UVC model, PHILIPS). By passing oxygen-feed gas through an ozone generator (O180F/DST, Canada), a steady flow of ozone gas was introduced into the reaction media by a micro-diffuser. O 3 concentration adjusting via flow rate of feed oxygen. According to 2350 E method (iodometric) of the standard methods for the examination of water and wastewater, ozone generator capacity and ozone concentration in the reactor and output were determined 26 . The contents of the reactor during the PCO/O 3 process was mixed by a magnetic stirrer. To measure OTC, 10 ml of each reaction mixture were filtered with 0.22 m syringe filters after being taken from the reactor and injected into high-performance liquid chromatography (HPLC). As the eventual results, mean values were calculated from three repetitions of each test. OTC, COD, and TOC degradation percentages in the PCO/O 3 process were calculated using the following (Eq. 16) 24 .
where, C o and C t are the concentrations of initial and residual OTC, COD, and TOC (mg/l) at time t, respectively.

Optimization of PCO/O 3 variables based on central composite design.
The undertaking of PCO/ O 3 process for OTC, COD, and TOC removal was optimized at 10 mg/l OTC as an initial concentration via central composite design (CCD) as an approach in the response surface method. The range of variables is illustrated in Table 1. Table 2 shows the number of 30 runs intended for the removal of OTC, COD, and TOC based on the CCD design.
The empirical second-order polynomial model of OTC removal via PCO/O 3 was described by Eq. (17) 27 :  www.nature.com/scientificreports/ Here, Y is OTC decrease (%),β 0 is the intercept, β i , β ii , and β ij are the linear, quadratic, and interaction effect coefficients of variables respectively, x i and x j are coded testing classes of the variables, k is the number of the independent variables, and e is the residual error. A coded value for each variable was calculated using Eq. (18), to enable comparison between factors with different units.
where, the coded value of the variable is Xi, ∆x which is the difference between the high and low values of the variable, X 0 represents the low value of the variable. Using P-values and F-values, analysis of variance (ANOVA) was performed to determine the interaction between response and factors. R 2 and R 2 adj are correlation coefficients. In addition, R 2 predicts were used. Then, In the optimal condition, organic, and inorganic radical scavengers effect, degradation kinetic based in difference OTC concentration, the synergistic effect on individual and combined processes (photolysis, simple ozonation, adsorption in dark mood, UV-O 3, photocatalysis, and PCO/O 3 process, the light source type (UV-C, UV-A, and visible light) effect, degradation pathway, and intermediates production, changes in wavelength scanning in during the PCO/O 3 process, and calculation of energy consumption were performed.
Analytical methods. Each sample was collected after the runs and processed through a 0.45 m microfilter before injection into HPLC. The amount of OTC was measured using a HPLC (Ciciel company, www.nature.com/scientificreports/ pump:120 mm × 240 mm × 440 mm, pressure: 40 MPa, flow rate setting range: 0.001-20 ml/min, solvent delivery method: parallel-type double plunger, plunger capacity: 10 μl) equipped with a UV-Vis detector (C18 column, 250 mm × 4.6 mm, with 5 μm particle size, pore size: 12 nm, surface area: 410 m 2 /g, carbon loading: 20%, pore volume: 1.25 ml/g, pH range 2-7.5, bonding type: monomeric) and a column. Temperatures of 4-25 °C and 353 nm, respectively, were kept constant for UV detection and column hold. Mobile phase consisted of 95:5 (v/v) mixture of acetonitrile and ultra-pure water (0.05% ammonia) at 1.0 ml/min flow rate. A 10 μl of OTC solution was injected at a retention time of 12 min. The TOC analysis was performed with Analytikjena's multi C/N 3100 TOC analyzer. As described in the standard method, COD was measured using the titrimetric approach (5220-C; closed-reflux) 28 . Residual catalyst and ROS were two interferers in COD test. So, the COD samples were initially passed through a 0.22 µ PTFE filter. Due to the rapid rate of decay, ROS are rapidly eliminated. The decay of the intermediates was studied by gas chromatography/mass spectrometry (GC-MS) (Agilent 7890A, California, USA). For the gas chromatography system, DB-5MS columns (30 m 0.5 lm film thickness) were used with highpurity (99.99%) helium as the carrier gas flowing at 1.0 ml/min. For temperature, the column was maintained at 35 °C for 1 min, increased to 300 °C at 7.0 °C/min, and held at 300 °C for 1 min. 10 cc of the sample was injected, with the injector and detector both set at 280 °C.
Supplementary studies. Kinetics of reaction and synergist effect. Degradation of OTC in the optimum condition of PCO/O 3 process described by the first-order reaction (Eq. 19) 29 .
where, C t is residual and C 0 is initial OTC concentrations (mg/l), t is the reaction time (min) and k app is the rate constant (1/min). The synergist effect of single mechanism such as photolysis, simple ozonation, catalytic ozonation, photocatalytic, and adsorption were investigated in an optimum condition of parameters. The synergist effect calculated via (Eq. 20) 30 . where P (kWh) is the consumption of power and V (l) is the volume of the treated solution.
Mineralization, by products, intermediates, and degradation mechanism studies. In this part of the study, we investigated mineralization, by products, intermediates, and degradation mechanism of OTC in optimum condition of PCO/O 3 process.

Results and discussion
Characters of MOF. XRD. The crystalline structure of the prepared BiOI, NH 2 -MIL125, and BiOI-MOF were analyzed through powder XRD patterns and the obtained finding are showed in (Fig. 2). The four major peaks in BiOI's XRD pattern (Fig. 2a) are (102), (110), (200), and (212) respectively. A tetragonal phase was observed in all diffraction peaks for synthesized BiOI (JCPDS no. 10-0445). High crystallinity is clear from the strong and sharp peaks. According to Han and coworkers' study, XRD analysis was similar 32 . MOF crystal formation is shown in (Fig. 2b). In the standard pattern, MOF peaks were observed at 6.7°, 9.8°, 11.9°, 15.06°, 16.6°, 17.9°, 19.07°, and 19.6°. It is clear from XRD patterns that orthorhombic phase correlates well with lattice parameters a = 15.03 Å, b = 6.2 Å, and c = 19.15 Å. By BDC and NH 2 BDC as a linker, three-dimensional pores are produced in a chine of corner-sharing TiO 6 . A similar result was reported by Ovisi and coworkers 33 . As shown in Fig. 2c, the BiOI-MOF structure also exhibits some peaks in the BiOI and MOF structures. According to the MOF, the prominent peak lies two degrees below 10 degrees. It appears at 25-35 degrees for BiOI peaks. Peaks' positions and intensity both change significantly. XRD results showed similarity in Du and coworkers' study (36). Based on Scherrer's equation (Eq. 13), the crystal size of BiOI, MOF, BiOI-MOF samples are 26/5, 47/8 and 37/1 nm.
FT-IR. In order to investigate the organic groups of samples, the FTIR spectrum was conducted (Fig. 3).
Five peaks at 571.7, 1309.5, 1382.2, 1627.4, and 3529.1 cm −1 position identified in BiOI (Fig. 3a). A 2 u type Bi-O bond has a symmetrical vibration at 571.7 cm −1 . BiOI exhibits strong adsorption at 1300 cm −1 to 1700 cm −1 , as well as a large peak at 3175 cm −1 to 3470 cm −1 as a result of flexural (-OH) and tensile vibrations (-OH) of free water molecules on the surface of the BiOI, respectively. Similar results reported by Ai and coworkers 34 . The peaks between 500 and 1500 cm −1 for MOF shows characteristics of organic compounds (Fig. 3b). Typical tensile Morphology analysis. Samples are analyzed by FESEM, EDX, EDS mapping and TEM (Fig. 4). FESEM images of BiOI, MOF, and BIOI/MOF was illustrated in (Fig. 4a). BiOI FESEM image shows the typical morphology that is thickness plate structures. It is relatively long to wide for a synthesized structure. An alternative solvent was used in Arumugam and coworkers' study. A 2D square-like nano-sheet morphology is observed with water (H 2 O-BiOI) as a solvent 37 . Synthesis methods, chemicals, and solvents used affect the morphology of BiOI. The synthesized MOF is thin and disk-like. A MOF's structure changes when its organic-metal ratio changes.
Research has shown that MOF structures change from rectangular to round as organic linkers increase. The nanoplate morphology of NH 2 -MIL125(Ti) was uniform in Zhang and coworkers' study 38 . BiOI-MOF is shown as a rod containing fine particles. The morphology of BiOI/MOF was similar in Du and coworkers' study 17 . It can be seen from the TEM image (Fig. 4c) that the NH 2 -MIL125(Ti) composite has been distributed successfully on the surface of BiOI, forming a core-shell structure. In the (Fig. 4c), shows the results of EDS-mapping. The results show that the distribution of elements on the surface of all three structures is heterogeneous. Also, (Fig. 4d) shows EDX image of BiOI/MOF. Samples show different weight ratios for each structure's fundamental elements.
BET analysis. BiOI/MOF, nitrogen adsorption-desorption isotherms are shown in (Fig. 5). Relative pressure affects adsorption and desorption rate. In view of the curved shape, the hysteresis belongs to category H 3 and    www.nature.com/scientificreports/ type IV. A non-hard cavity with an incised shape has this type of hysteresis. The tensile strength effect, also known as the H 3 hysteresis repulsion branch, has a steep slope. By N 2 adsorption-desorption isotherm and BJH, the precursor's BET area of surface and pore volume were determined to be 947.85 m 2 g and 16.27 cm 3 /g, respectively.
UV-Vis analysis. Samples UV-Vis spectrum shown in (Fig. 6). The band edge of the BiOI, MOF and BiOI-MOF calculated at 635, 505 and 685 nm, respectively. The results show that the catalyst BiOI/MOF has a larger band edge compared to pure BiOI and MOF. The band edge is in the visible light range, indicating that a visible light source can also activate the catalyst. Also, DRS analysis showed that the BiOI/MOF catalyst has a shorter band gap than BiOI and MOF and requires less energy for activation. In our pervious study, similar results were reported 26 .
XPS analysis. The chemical state and surface composition of BiOI-MOF were determined using XPS (Fig. 7). The XPS spectrum confirms elements O, Ti, C, I, and Bi. The composition purity was verified by BiOI and NH 2 -MIL125(Ti). An overview of BiOI-MOF structure elements is presented in Table 3. Similar results were obtained by Jiang and coworkers 18 .
Photocatalytic-ozonation activity. CCD design, fit, and statistical analysis. The statistical analysis of the OTC, and COD results provided the basis for response surface method (RSM) development and finding the optimal conditions for PCO/O 3 process ( Table 2). F-value and P-value in the analysis of variance (ANOVA) analysis are substantiated statistical indices for adjusting how to fit the data deviation factors. Because of these indexes, a significant statistical model with a high F-value and a low p-value (0.05), was selected. For the relationship between the predicted and experimental values of OTC and COD, Fisher's F-test indicated that a quadratic model would be the most appropriate fit (Table 4). www.nature.com/scientificreports/ Using the ANOVA, the results of OTC and COD (Table 5) appeared highly reliable and had a very low probability value for the quadratic regression model, suggesting that it could accurately explain the codes within the actual data and predicted values. Based on the correlation coefficients (R 2 , R 2 adj, and R 2 predict), this model proved exceptionally valid for predicting responses. There is a perfect correlation between R 2 , R 2 adj, and R 2 predict, which suggests that the model has a high ability to predict responses 39,40 .
The F-value of OTC, and COD is 2451/04, and 450/34 respectively. In the OTC and COD model, P-values of 0.0001 indicate a robust fit to the relationship between the experimental and predicted response values. A model discernment greater than 4.0 for the Adequate Precision ratio indicates a satisfactory model 41     www.nature.com/scientificreports/ A positive and a negative effect of the studied mediators on the removal of OTC and COD is illustrated in Eqs. (22) and (23). The quadratic model had a low coefficient of variation (CV < 10%), which confirms the precision and reproducibility of the experiments. All interaction influences among all variables were significant. Based on the statistical analysis, the correlation between variables and the removal of OTC and COD was found to be statistically significant. (23) COD removal (%) = 72 + 6/3 × A + 4/1 × B + 7/5 × C + 1/2 × D − 0/06 × AB + 0/18 × AC + 0/4 × AD www.nature.com/scientificreports/ The impact of parameters on the OTC and COD removal on PCO/O 3 process. PCO/O 3 process highly depends on the pH of the solution, which has a substantial impact on the removal efficiency. Solution pH affects pollutant hydrolysis, ionization pollutants, catalyst surface characteristics, and activity rate of oxidants and reactive species, as well as the degradation route. Hence, the effects of pH (4-8) on OTC and COD removal efficiency were examined and illustrated in (Fig. 8). OTC and COD removal efficiency significantly increased with increasing solution pH, with pH of 8.0 achieving the highest removal efficiency. Electrostatic interactions between catalyst functional groups can be affected by pH (e.g., ionization) in this case. Also, when pH increases, ozone decomposes rapidly into radicals, which can degrade organic compounds much more rapidly and efficiently than ozone molecules 42 . It has been reported that OTC molecules have two pka values (3.18 as the strongest acidic and 8.29 as the strongest basic); therefore, OTC molecules are in cationic and anionic forms at the pka below 3.18 and above 8.29, respectively 43 . Meanwhile, the pH zpc of the photocatalyst was 6.3, indicated that the photocatalyst surfaces act as protonates and non-protons in lower and higher, respectively. At lower pH levels, minimal adsorption was observed because of electrostatic repulsion between the positively charged catalyst and the protonated OTC molecules. Many studies have shown that contaminant decomposition occurs primarily through direct and indirect oxidation pathways in the ozonation process. OTC and ozone react directly at lower pH levels. As the pH of the solution was increased, ozone decomposition rate and free radical production increased, and free radicals with high oxidation potential were generated. These radicals demonstrate superior ability in OTC degradation. OTC molecules can react indirectly with free radicals in an alkaline environment, thus resulting in a higher removal efficiency. Asgari et al. investigated the ability of the photocatalytic ozonation process for ceftazidime removal. Based on this research, pH value at 11.0 was the optimum pH of solution 44 . The similar findings have been also reported by Espíndola et al. 45 . In this study, pH equal to 7.5 was reported as an optimum pH solution to oxytetracycline oxidation in the photocatalysis process.
Dose of photocatalyst. The dose of the photocatalyst is another important parameter that affected heterogeneous AOPs based on the catalyst. Different amounts of BiOI/MOF were applied to determine PCO/O 3 efficiency in OTC removal (Fig. 8). Based on results, the removal efficiency of OTC was significantly enhanced by increasing the catalyst dose from 0.25 to 0.5 mg/l over 10 ppm of the initial concentration of OTC. Nevertheless, 0.4 mg/l catalyst doses did not show any notable differences. Higher catalyst concentrations had less favorable effects on the PCO/O 3 process. Typically, this phenomenon is caused by particles accumulating and agglomerating and, reducing their surface active sites. Excessive catalyst dosages can scatter UV light in the solution, preventing the transmission of ultraviolet light and ozone to the catalyst surface. The surface area of the active site for photocatalytic activity can be increased by increasing catalyst content. As a result, a greater amount of ozone was adsorption on the surface of the nanoparticles, resulting in the production of more active radicals as a result of the breakdown of ozone. This improved OTC removal efficiency. Nanoparticles produce oxygen radicals when combined with ozone under alkaline conditions. When oxygen radicals are present in water, they create · OH radicals, which eventually enhance ozonation efficiency 42 . Photocatalytic ozonation studies have reported similar findings, which agree with this study. Here, Lu et al. study tetracycline hydrochloride degradation by photocatalytic ozonation. In this research, zero to 0.8 g/l Bi 2 WO 6 was applied as a catalyst. Until 0.5 g/l, performing process increased to 78% in 120 min then at 0.8 g/l the efficiency dropped to 65% 46 . Also in Yu et al. study, similar result obtained 47 .
Ozone concentration. Main important economic parameter in PCO/O 3 process is ozone concentration. This parameter is closely related to energy consumption. An ozone gas concentration of 20-40 mM/l-min was investigated in this study (Fig. 8). The results showed that by increasing the concentration of ozone gas from 20 to 30 mN, the process's efficiency increased and decreased slightly. Higher ozone concentrations increase mass transfer in the reaction media. Concentrations of dissolved ozone in the solution increase as the ozone gas flow rate increases. Because of a synergistic effect, reactive oxygen species (ROS), particularly hydroxyl radicals, are produced in greater quantities. Photocatalyst easily adsorb dissolved ozone molecules because of weak hydrogen bonds with their surface · OH groups. Anions of ozonide radicals are produced by capturing electrons on the surface of the photocatalyst. Thus, the ozonide radicals (O 3

·−
) increase the amount of · OH radicals and the efficacy of OTC degradation 46 . When incoming air flow increases excessively, and the mass transfer rate from the gas phase to the liquid phase is limited, the efficacy of removing decreases. The concentration of ozone required for PCO/ O 3 depends on the type of process, the type of reaction reactor, the type of pollutant, and the specifications of the intermediate compounds 48 . In previous studies, a different concentration of ozone was reported as an optimum concentration. Yu et al. 49 1.5 mg/l-min of ozone concentration, Heydari et al. 50 11 mg/l of ozone concentration was reported as an optimum ozone concentration.
Reaction time. The reaction time is another important parameter that affected heterogeneous AOPs. In the current study, the effect of reaction rime on PCO/O 3 effect in OTC and COD removal in 30 to 60 min range was investigated (Fig. 8). Based on (Fig. 0d), significant increase in degradation efficiency is observed for both OTC and COD with increasing reaction time. As mentioned, PCO/O 3 process removes OTC through two direct and indirect oxidation mechanisms. Direct oxidation by ozone molecules and photolysis and indirect oxidation based on free radicals and electron-hole pairs produced cause pollutant destruction. So increasing the reaction time allowed for complete pollutant destruction and mineralization, and less intermediate compounds formed. Therefore, the optimal reaction time depends on the type of pollutant, characteristics of the reaction reactor, other operating parameters, and operating conditions. For this reason, different times have been reported in previous studies. Lu et al. 46 , report 120 min as the optimal time. Synergist effect, TOC removal and recyclability of photocatalyst. Figure 9a shows the results of comparing the roles of the single, binary, and PCO/O 3 processes in the degradation of OTC under optimal conditions. According to the results, simple ozonation (55%), photolysis (UV-C = 30%, UV-A = 12, and visible light = 8%), and sorption (7%) are ineffective in OTC degradation. Therefore, these mechanisms do not have satisfactory potential for degrading OTC individually. O 3 /catalyst (68%), UV/ catalyst (60%), and UV/O 3 (66%) as binary processes exhibited more efficient OTC degradation than the individual processes. Ultimately, over 95% of the tetracycline concentration was removed through the fundamental process. This study confirms a synergistic effect between various processes, resulting in a more active degradation of OTC and an excess creation of super reactive radicals. However, when all variables were involved simultaneously, OTC removal performance significantly improved. Based on (Eq. 20), the SF coefficient of 1.31 was obtained. When light irradiates to catalyst, electrons are produced in the conduct band (CB) and holes are formed in the valance band (VB). Since the catalyst's CB potential is higher than the redox potential of oxygen/oxygen peroxide in the reaction medium, electrons transfer to the oxygen in the reaction media and are converted into oxygen abundantly. The photocatalytic system will be affected by the utilization of the photo-generated carriers and the production of ROS. In contrast to its direct ozonation, electrophilic ozone easily traps photo-generated electrons when introduced into photocatalytic systems. Because of the effective transfer of the photo-generated electron of catalyst, the electron-hole pairs www.nature.com/scientificreports/ are efficiently separated under simulated light irradiation, leading to a large amount of ROS being produced. A catalyst can also improve the O 3 /Light system's ozone utilization ratio. In the reaction, more ozone will produce more ROS, which is beneficial to mineralization in removing organic compounds. It is possible to remove organic pollutants effectively by combining photocatalysis with ozonation, as the combined process produces more ROS than the sum of the individual processes. In Lu et al. 46 study, TOC removal in O 3 , O 3 /light and O 3 /light/catalyst were 18, 20 and 78% respectively. In (Fig. 9b) the results of OTC, COD, and TOC removal in 0-100 min reaction time presented. According to the results, it has been found that the COD and TOC removal rate is lower than OTC. Under optimal conditions, OTC, COD, and TOC decreased by 96.2, 77.2, and 64.2%, respectively. The existing difference indicates the formation of organic intermediates and by-products. Also, the results showed that increasing the reaction time increases the efficiency of removing COD and TOC and more mineralization. Asgari et al. 42 reported 100, 92.3 and 81.1% CIP, COD and TOC removal respectively, in 60 min. Stability and recyclability are essential characteristics of catalysts for real applications. As a result, photocatalyst recyclability was investigated over six consecutive cycles under optimal conditions. After each cycle, the catalyst was separated from the solution and dried for 24 h at 80 °C. According to (Fig. 10 , OH · , e --h + ) were evaluated by quenching tests for their effect on removing OTC 52 . OA acts as a scavenger of O 2 ·and e --h + pairs. TBA is the scavenger of O 2 1 and OH · , and finally N 2 gas is a scavenger of all ROS s . Based on results indicated that, performing PCO/O 3 process in the presence of AO, it drops less, in the presence of TBA, the rate of decrease in process efficiency is higher, and N 2 gas causes a significant decrease in efficiency. Therefore, O 2 1 and OH · were the main active radicals. Similar finding reported in previous studies 48,53 .
Kinetic of reaction and electrical energy consumption. On optimized conditions, the pseudo-first-order model (Eq. 19) described the OTC degradation rate over the PCO/O 3 process. As shown in (Fig. 12), OTC degrades over PCO/O 3 in 10, 20, and 30 mg/l concentrations over a 60 min reaction using a linear pseudo-first-order kinetic model. According to the pseudo-first-order model, OTC degradation kinetics in PCO/O3 are well supported by the experimental data (R2 > 0.99). Based on the review of previous studies, it was found that in most AOPs, the reaction kinetics follows pseudo-first-order kinetics 54,55 . www.nature.com/scientificreports/ In practical applications, energy consumption must be considered. The ozone generator and the UV lamp handled EEC in this process. The total energy consumption of was 220 W. According to (Eq. 21), EEC of PCO/ O 3 process under optimum condition was 4.44 KWh. Based on previous studies, the EEC depends on various factors such as the power of the ozone generator, the power of the lamp, the reaction time, the volume of the reaction chamber and finally the efficiency of the process 56 . In Kang et al. study, 55 KWh/m 3 reported as EEC 10 .
Identifying the mineralization pathway for OTC. To identify the intermediates produced during OTC degradation, GC-MS experiments were conducted to reveal the degradation pathway. According to the results analysis and previous relative studies, two main pathways were identified for OTC degradation: direct oxidation and indirect oxidation. During direct oxidation, mineralization plays a minimal role. Indirect oxidation destroys OTC through hydroxylation, demethylation, and decarbonylation (Fig. 13). OTC's aromatic ring is attacked by hydroxyl radicals (pathway 1) through the reaction of the OH addition. By-product P1 had been formed after OH attacked the keto-Enol at C11a-C12 (m/z 476). As a result of the decarbonylation (pathway 2), OTC  www.nature.com/scientificreports/ loses CO from its ring structure, which produces the byproduct P6 (m/z 433). In the chemical decomposition of acetylacetone, CO was generated by the rapid dissociation of COCH 3 formed by the cyclic cleavage of enolic acetylacetone. At the C4 target site, one methyl group of dimethylammonium group might be eliminated, resulting in intermediate P3 (m/z 448) 57 . By review of previous studies, it was found that based on the process used, different intermediate compounds and different reaction paths may occur in the degradation of OTC.
To conclude, Table 6 demonstrates multiple studies that assessed different AOPs strategies for removing OTCs from aquatic solutions.

Conclusion
In this study, BiOI-MOF was synthesized by solvo-thermal method. Based on results of XRD, FESEM, EDAX, FTIR, UV-Vis spectra, TEM, BET, and XPS analysis indicated that synthesis of BiOI/MOF was done with best character. Design of experiment, ANOVA statistical analysis, interaction of parameters and predicated optimum condition was done based on CCD. The effect of pH of solution (4)(5)(6)(7)(8), catalyst dose (0.25-0.5 mg/l), reaction time (30-60 min) and O 3 concentration (20-40 mN) at 10 mg/l of OTC on PCO/O 3 process optimized. Based on F-value and P-value coefficients, the model of OTC (F-value = 2451.04, P-value = 0.0001) and COD (F-value = 450.3, P-value = 0.0001) removal was quadratic. Under optimum condition pH 8.0, CD = 0.34 mg/l, RT = 56 min and O3 concentration = 28.7 mN, 96.2 and 77.2% of OTC and COD removed, respectively. TOC decreasing under optimum condition was 64.2% that lower than OTC and COD. The kinetic of reaction follows pseudo-first-order kinetics (R 2 = 0.99). Synergistic effect coefficient was 1.31 that indicated ozonation, presence of catalyst and photolysis had a synergistic effect on OTC removal. The stability and reusability of catalyst in six consecutive operating was acceptable and 7% efficiency decreased only. Cations (Mg, and Ca), SO 4 2− had no effect on the efficiency of the process, but other anions, organic scavengers and N 2 gas, had an inhibitory effect on the efficiency of the process. EEC under optimum condition 4.44 KWh/m 3 calculated. Finally, the OTC degradation probably pathway includes direct and indirect oxidation that Hydroxylation, Demethylation, and Decarboxylation were the fundamental mechanism in OTC degradation.

Data availability
The datasets generated and analyzed during the current study were available from the corresponding author on reasonable request.